Probing the excitation spectrum of nonresonantly pumped polariton condensates 
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A series of remarkable experiments have recently 
demonstrated the occurrence of Bose-Einstein conden- 
sation in systems of exciton-polaritons in semiconductor 
microcavities [H, [3, H, 0, HI • These observations, together 
with parallel ones on magnon condensation in magnetic 
solids (6|] , are opening the way to the investigation of the 
Bose-Einstein condensation phase transition and of Bose- 
condensed quantum fluids in completely new regimes far 
from thermal equilibrium where the state of the conden- 
sate is no longer determined by a thermal equilibrium 
condition, but by a dynamical balance between driving 
and dissipation [7]. 

A key quantity in the theory of quantum fluids is the 
dispersion of the elementary excitations, which deter- 
mines the dynamics of the system in response to external 
perturbations and, in particular, plays a central role in 
determining its superfluidity properties [8j]. While a de- 
tailed knowledge is nowadays available of the phonon and 
roton branches of liquid Helium and of the Bogoliubov 
modes of dilute atomic gases [9|, not much experimental 
work has been performed yet on polariton condensates. 
On one hand, pioneering luminescence [l^ and pump- 
and-probe experiments [ll|, |l2] have addressed the ex- 
citations of resonantly pumped condensates in an OPO 
configuration; on the other hand, luminescence experi- 
ments suggesting a linear dispersion of the elementary 
excitations in nonresonantly excited polariton conden- 
sates have been recently reported [l^. Still, none of these 
works has provided complete evidence of the peculiar fea- 
tures that were predicted to appear because of the non- 
equilibrium nature of polariton condensates [1^, ll5|, ll6| . 

From a different standpoint, while quite some evi- 
dence is available for "off-branch" scattering processes 
in OPO configurations [lO^Ilii, no observation of the re- 
lated negative energy "ghost" branch that is expected to 
appear under non-resonant pumping as a consequence of 
polariton-polariton interactions has been reported yet. 
Most likely, this stems from the fact that exp eriments 
were based on luminescence spectroscopy [13|, so that 
the emission from the ghost branch was easily masked by 
the much stronger background of the condensate emis- 
sion. A few recent theoretical work has proposed more 
refined schemes that appear suited to overcome this dif- 



ficulty and detect this elusive branch by looking at either 
the absorption [17] and resonant Rayleigh scattering fl^ 
spectra, or at the density response to an external pertur- 
bation [19]. 

In the present work, we push these ideas forward and 
we propose a simple four-wave mixing (FWM) spec- 
troscopy scheme to measure the dispersion of all elemen- 
tary excitations branches of a non-resonantly pumped po- 
lariton condensate: of the three incident beams needed in 
a FWM experiment, two are provided by the condensate, 
and one by an external laser field. The dispersion of ex- 
citations is inferred by scanning the energy and wavevec- 
tor of this latter beam and recording the transmitted, 
reflected, and/or four- wave mixed beams. A first appli- 
cation of related FWM techniques to polariton systems 
across the parametric oscillation threshold was reported 
in [2^ , but no specific interest was paid to the excitation 
modes of the condensate itself. 

The mean-field Gross-Pitaevskii formalism to describe 
the polariton condensate dynamics that was introduced 
in [16] is briefly reviewed in SecU and then applied in 
Secinito describe the FWM response of polariton conden- 
sates. Closed formulas are extracted for the transmission, 
reflection and FWM signals in the spatially homogeneous 
case and then used to discuss the main features of the 
spectra. Generalization to the experimentally relevant 
case of finite-size condensate is investigated in SeclIIII 
the efficiency of the FWM method for realistic conden- 
sates is confirmed by means of a numerical solution of the 
polariton GPE in non-uniform geometries. Conclusions 
are finally drawn in Sec. [iVl 



I. THE ELEMENTARY EXCITATION 
SPECTRUM 

As a consequence of the short life time of polaritons in 
state of the art microcavities, the polariton condensate 
can hardly be considered as a thermal equilibrium ob- 
ject: Continuous external pumping is necessary to keep 
the condensate in a stationary state, and the state of 
the condensate is determined by a dynamical balance 
between pumping and dissipation. While the momen- 
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turn distribution of large systems still appear to follow 
at large energies the typical exponential law of equilib- 
rium statistical mechanics [H, [3, 0] , clear evidence of the 
non-equilibrium nature of the polariton condensate has 
been observed in the ballistic outward polariton flow from 
small-sized condensates [2l|, [22| as well as in the sponta- 
neous appearance of vortices in the presence of a signifi- 
cant disorder potential [23]. 

Several theoretical papers have recently investigated 
the effect of the non-equilibrium condition on the elemen- 
tary excitation spectrum [13, [IB] and in particular have 
pointed out the diffusive nature of the Goldstone mode 
at low wavevectors. In what follows, we shall adopt the 
point of view of our previous paper ^] where the dynam- 
ics of a non-equilibrium condensate was discussed at the 
mean- field level in terms of a generalized non-equilibrium 
Gross-Pitaevskii equation. 




FIG. 1: Plot of the elementary excitation spectrum in a 
spatially homogeneous geometry and sketch of the proposed 
FWM experiment. The red dot indicates the wavevector 
A; = and the energy uoc of the condensate. The solid green 
line is the dispersion ([T]) of the elementary excitations on top 
of the condensate. The arrows show the parametric scattering 
process on which our proposed scheme is based: the conden- 
sate is probed at the lower green dot and the response is 
probed via the emission at the upper green dot. The dotted 
black line is the dispersion of free polaritons at linear regime. 
The solid and dashed red lines include the blue-shift due to 
condensate-condensate and condensate-reservoir interactions 
as mentioned in the text. System parameters: 7 = F = gric^ 
gnriR = ^.bgric- 

This gives the following form for the dispersion of the 
elementary excitations on top of a spatially homogeneous 
non-equilibrium condensate: 

uj^{k) -ujc = -iT/2 ± ^[c^Bo,(fe)]'-rV4 (1) 

which is plotted in FigH] as a green solid line. In the 
low-/c, diffusive region of the spectrum the ± signs corre- 
spond to respectively the phase and the density branches. 



In the high-/c region where the standard equilibrium Bo- 
goliubov dispersion is recovered, they instead correspond 
to the normal and anomalous "ghost" branches. In the 
following of the paper, a specific attention will be paid 
to the "ghost" branch which is a signature of coherent 
polariton-polariton interactions. 

The condensate energy uOc (red dot) is blue-shifted 
with respect to the bottom of the lower polariton branch 
(black dotted line) by the condensate-condensate and 
condensate-reservoir interactions, c^c = ^0 + + 9r^r- 
Here, ric and ur are the densities respectively of the con- 
densate and of the exciton reservoir, while the g and 
gR coupling constants characterize the interactions re- 
spectively between a pair of condensate polaritons and 
between condensate polaritons and reservoir excitons. 
The Bogoliubov spectrum of equilibrium condensates 
has the standard form ujBog{k) = [ski^k + in 
terms of the bare polariton dispersion at linear regime. 
Within the parabolic approximation, this can be written 
as s/c = hk'^ /2mLP' As discussed in the effective 
width F is defined asF = a/37/(l + a/^), in terms of the 
scaled pump intensity a = P/ Pth — 1 above the threshold 
value Pth and the 13 coefficient defined in terms of gain 
rate R{nR) of the condensate from the excitonic reservoir 
SiS (3 = n%R' {n%) /R{n%). Well above the threshold 
a 1, the effective width F recovers the empty-cavity 
polariton decay rate 7. 

Note how the dispersion of elementary excitations 
significantly differs from the naive Hartree prediction 
CO = coq -\- 2gRnR + gric + hk'^ /2m (dashed red line) and 
rather recovers in the high-momentum region the value 
uj = uJo-\- gn^R + 2gnc + hk'^/2m which correctly includes 
a factor 2 due to bosonic exchange (solid red line). Pro- 
vided the effect of saturation [25| on the polariton inter- 
actions is negligible, these simple arguments can be used 
to isolate the contribution of respectively the condensate- 
reservoir and the condensate-condensate interactions to 
the blue-shift, and possibly even evaluate the fraction of 
condensed excitons with respect to the total number of 
excitons in the microcavity. 



II. THE FOUR WAVE MIXING 
SPECTROSCOPY SCHEME 

Even though a simple luminescence experiment should 
in principle be able to reveal both the positive and the 
negative "ghost" branch of the elementary excitation 
spectrum shown in Fig{Tl experimental observations have 
not been able to provide unambiguous evidence of the 
"ghost" branch yet [l|, 0, |l^ : most likely the very weak 
luminescence coming from the "ghost" branch is hidden 
by the much stronger emission from the condensate and 
the upper branch. 

The physical process underlying our proposal is 
sketched in Fig. [H Elementary excitations are created 
on top of the condensate by injecting extra polaritons 
with a probe laser beam at a finite in-plane wavevector 
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k tuned at a frequency uj. The response of the system 
is then observed via the coherent hght emission at an 
opposite wavevector — k at an energy 2uOc — oo: the exis- 
tence of a coherent coupling between the symmetrically 
located frequencies uo and 2uOc — oo and wavevectors ±k 
stems from the fact that the elementary excitations of the 
condensate consist of a coherent superposition of plane 
waves at (k, cj) and (— k, 2a;c — oo) [9, 26]. Equivalently, 
this same process can be interpreted in nonlinear opti- 
cal terms as a stimulated parametric scattering where 
a pair of condensate polaritons are coherently scattered 
into one more probe polariton plus one FWM polariton 
of symmetric wavevector and frequency. Within the stan- 
dard language of four- wave mixing, two of the three input 
beams are provided by the condensate itself, and only the 
third comes from the incident laser. 



A. Transmittivity 

The same linearization procedure that was used in [l^ 
to obtain the dispersion of the elementary excitations can 
be straighforwardly used to determine the response of the 
system to a (weak) probe laser of amplitude £^pr(k, cj). 
The simplest quantity to consider is the probe transmis- 
sion. 

In the usual approximation that the effective reservoir 
relaxation rate is much larger than all other frequency 
scales in the problem, the reservoir dynamics can be elim- 
inated and the amplitude of the transmitted field has the 
form: 



/To 7 gnc^e\^ 



(cj-c^+(k))(c^-c^_(k)) 



(k,c^). 



(2) 

The frequencies a;±(k) that appear in the resonant de- 
nominators were defined in ([T]): the transmission is res- 
onantly enhanced when the the probe is on resonance 
with an elementary excitation of the polariton conden- 
sate, either on the positive energy branch or on the neg- 
ative energy, ghost branch. A complete density plot of 
the transmittivity T = \Etr/Epr\^ in the (cj, /c) plane is 
shown in Fig|2ja). For the sake of simplicity, the trans- 
mission T is normalized to the linear regime resonant 
transmittivity Tq of the unloaded cavity. 

At large momenta /c, the intensity of the ghost mode 
is much weaker than the one of the upper, normal mode 
and quickly tends to zero. Already for /c^ = 1, the ghost 
resonance is extremely weak and almost invisible in the 
cut shown as a dotted line in FiglSl^c). On the other hand, 
the normal resonance tends to the unloaded cavity peak 
of height T/To = 1 and width 7. For decreasing values 
of /c, the intensities of the two resonances get closer and 
eventually stick at the singular point /c^ = 0.5 of the 
dispersion [Fig|2l^e)]. 

In the low-/c diffusive region where the two branches 
are degenerate Re[c(;±] = 0, the transmission intensity 
significantly exceeds the peak transmittivity Tq of the 
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FIG. 2: Transmittivity (a) and FWM (b) spectra of a non- 
equilibrium condensate excited with a probe laser at wave 
vector /c > and frequency cj. The full and dashed green lines 
indicate the positive and negative Bogoliubov branches from 
Eq. ([1]). Cuts of the transmission and FWM spectra at given 
k values are shown in panels (c) and (d); the corresponding k 
values are indicated by the vertical lines in panels (a,b). The 
peak value of the transittivity and the FWM signal on the 
respectively the normal and ghost Bogoliubov branches are 
shown in (e) and (f) as solid and dashed lines. In (f) the two 
quantities coincide. All spectra are normalized to the empty- 
cavity, linear regime transmittivity Tq. System parameters: 
F = 7 ^ gric- 



unloaded cavity [solid line in Fig|2fc)]: this gain feature 
is due to the condensate, that is able to amplify the probe 
beam. A similar phenomenology was found and discussed 
in ^] for the case of an optical parametric oscillator. 



While the appearance of the ghost branch in the trans- 
mission spectra of Fig|2] is a signature of the presence of 
the condensate, the transmittivity on the normal branch 
can exceed the linear regime value Tq even for pump in- 
tensities well below the condensation threshold. Even 
if not sufficient to overcome losses, some gain is in fact 
present also in this case [17], as witnessed by the polari- 
ton linewidth decreasing below the unloaded cavity value 
7- 
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B. FWM signal 

The same method can be used to obtain an expression 
for the amphtude of the four wave mixing signal: 



to7 



gric + iT/2 



{lj — a;+(k))(a; — a;_(k)) 



(3) 



It is interesting to note that both the collisional (oc guc) 
and gain saturation (oc F) nonlinearities contribute to 
the FWM signal. The response function SpwM = 
\Efwm / Epr\^ is plotted in the upper- left panel of Fig 121 
Thanks to the symmetry under the exchange uj —uj^ 
the heights of the two peaks corresponding to the normal 
and the ghost branches are equal. 

This is a crucial advantage of the FWM technique as 
compared to the transmission spectroscopy discussed in 
Sec lII Al or the resonant Rayleigh scattering discussed 
in Differently from this latter scheme, FWM does 

not rely on the presence of a disorder potential and there- 
fore is not affected by speckle-like modulations of the de- 
tected signal. As it happened for the probe transmission, 
also the FWM signal results strongly enhanced in the dif- 
fusive region and may eventually become stronger than 
the probe itself. On the other hand, it decreases quite 
quickly for larger values of k. 

Note that both direct polariton-polariton collisions and 
nonlinear gain saturation effects contribute to the FWM 
signal via respectively the terms g ric and ir/2 in the nu- 
merator of (|3]) : some FWM signal is therefore expected to 
appear even for negligible polariton-polariton collisions, 
but is in this case peaked very close to the free polariton 
dispersion. On the other hand, the presence of a FWM 
signal appears to be a conclusive evidence of the presence 
of a coherent condensate: the same calculation below 
threshold would in fact give a vanishing coherent FWM 
amplitude in the weak probe limit. Gain saturation is 
in fact not active in this regime and FWM processes in- 
volving the exciton reservoir only provide an incoherent 
background without any specific resonant feature. 



C. Reflectivity 



Note that the linewidth jph that appears at the numer- 
ator of (IH) is the linewidth of the cavity-photon mode in 
the absence of any excitonic resonance. This quantity is 
to be distinguished from the linewidth 7 of the polariton 
branch that appeared in the previous formulas. In terms 
of the Hopfield coefficients U^jf^ that quantify the cavity- 
photonic and the excitonic contents of the lower polariton 
branch, the polariton 7, exciton and cavity-photon jx,ph 
linewidths are related by 



LP\2 



(5) 
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FIG. 3: Panels (a,b): {k,oo) plot of the reflectivity for 7 = 
r = gric (a) and F = O.47 = gric (b). (c,d) show reflection 
spectra at fixed values of k as indicated by the vertical lines 
in panels (a) and (b). The peak values of the reflectivity on 
the normal and ghost Bogoliubov branches are respectively 
shown in (e) and (f) as solid and dashed lines. Parameters: 



As many microcavity samples are grown on an absorb- 
ing substrate, transmission measurements are not always 
possible. While the FWM expression Q is the same in 
both reflection and transmission geometries, the inter- 
ference between direct reflection of light on the external 
mirror and the secondary emission from the cavity makes 
the expression of the reflection amplitude a bit more com- 
plicate than the transmission amplitude (|2]): 

TP n ^ -Tp^ + Sk + ^ + >^ 

V 2 {uj - u;^{k)){u; - u-{k)) J 

(4) 



As it involves subtle interference effects, the reflection 
spectrum turns out to be sensitive to the ratio of the 
linewidth of excitations in the empty and pumped cavity 
7/r > 1, as well as to the ratio between the cavity-photon 
and the polariton linewidths jph/j- Here we focus our 
attention to the experimentally most relevant case where 
Iph = Ix = 1' Examples of spectra are shown in Fig. [3l 
for different values of the pumping intensity. Panels (a,c) 
refer to a case well above threshold for which r/7 1, 
while panels (b,d) are made for a lower (but still above 
threshold) value of the pump intensity at which r/7 = 
0.4. 



In panel (a), a minimum in the reflectivity is visible 
on resonance with the positive Bogoliubov branch, while 
the reflected intensity on resonance with the negative 
branch exceeds the incident intensity. As for the trans- 
mission and the FWM signal, the ampliflcation effect is 
the strongest in the diffusive region at low k [panel (c)]. 
At high values of /c, the dip corresponding to the normal 
branch remains fully visible, while the peak correspond- 
ing to the ghost branch disappears. This once again con- 
flrms that the FWM mixing scheme is the most suited 
tool for the detection of the ghost branch. 

The same quantities are plotted in panels (b,d,f) for 
a smaller value of r/7 = 0.4, i.e. for a pump intensity 
closer to the threshold. As a consequence of the narrower 
resonance line, the reflected beam is now more intense 
than the incident one on both the normal and the ghost 
branches. The ampliflcation in the \ow-k diffusive region 
is also stronger than in the case P ^ Pth considered 
before. 



III. PROBING THE ELEMENTARY 
EXCITATION SPECTRUM IN FINITE SIZE 
AND/OR DISORDERED CONDENSATES 

Even in the absence of trapping potentials, the spa- 
tial extension of polariton condensates is generally lim- 
ited by the size of the laser spot that is used to pump 
the microcavity 0. Typically, state of the art polari- 
ton condensates have a typical size of the order of tens 
of microns and may show a considerabe inhomogeneous 
broadening of their spectral features as a consequence of 
the inhomogeneous density proflle. Although this effect 
can be reduced by choosing top-hat pump beams, the in- 
terpretation of FWM experiments is the clearest if one 
focusses the probe spot onto a small region at the center 
of the condensate where the density is almost flat and 
one spatially selects the emission from the same central 
region. 

The result of numerical calculations based on the full 
polariton Gross-Pitaevskii equation is summarized 
in FigHl The case of a probe beam of diameter 1 jiva 
much smaller than the full condensate size and tuned at 
a frequency slightly below the condensate frequency uOc 
is considered. 

The case of a quite high probe wavevector kp = 
1.5/im~^ tuned close to resonance with the ghost branch 
is shown in the right panel of the flgure. As the k distri- 
bution is narrow as compared to the central wavevector 
kp , the response closely resembles the one of the spatially 
homogeneous case discussed in the previous section. A 
single FWM response peak appears around —kp and with 
a /c-space linewidth narrower than the one of the incident 
probe: only those wavevector components that are on 
resonance with the ghost branch effectively contribute to 
the FWM signal. 

A remarkable feature of the ghost branch is visible in 
the real space density proflle of the probed condensate 
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FIG. 4: (Color online) Numerical results for the (/c, uj) emis- 
sion pattern from a probed polariton condensate in a trans- 
mission geometry. The condensate frequency is Uc = 0.83 
meV and probe laser has a l/xm diameter. The central wave 
vector and the frequency of the probe are kp = 0.6 /xm"^, 
iUp — 0.3 meV (left panel) and kp — 1.5 /xm"^, cOp — —0.2 
meV (right panel) . The full and dashed lines show the analyt- 
ical prediction the excitation spectrum and the linear regime 
polariton dispersion. The central surface plot shows the real 
space polariton density profile corresponding to the right hand 
panel. In all panels, the condensate is excited with a top hat 
laser with a 30 /xm diameter and only the emission originat- 
ing from a central region with a 20 /xm diameter is recorded 
in the left and right panels. The finite frequency spread of 
the emission lines originates from the finite time window of 
about 12 ps from which the spectrum was extracted. 



shown in the central panel of the flgure: as the group 
velocity of the ghost branch is opposite to the wavevec- 
tor k, the perturbation due to the probe laser (with a 
central wavevector kp > pointing to the right) propa- 
gates through the condensate in the leftwards direction 
and concentrates on the left-hand side of the probe spot. 

The case of a probe laser with a smaller wave vector 
and a frequency closer to, but still below the condensate 
frequency is shown in the left panel. This choice dra- 
matically modifles the qualitative shape of the response: 
in particular, note how the maximum of the response 
at the probe frequency lies between the two excitation 
branches of the homogeneous system, a spectral region 
where collective modes are strongly affected by the flnite 
size of the system and possibly by the non-trivial spectral 
shape of the condensate [22^, [27] . Although information 
on the excitation modes of the condensate is hardly ob- 
tained in this conflguration, the very presence of a FWM 
signal is still a direct proof of the presence of a coherent 
condensate. 

The situation is a bit more complex in the presence 
of disorder. As any experimental sample is inevitably 
far from being perfectly homogeneous, it is important 
to briefly summarize the main features that we observe 
in our numerical simulations when a disorder potential 
acting on polaritons is included. A complete study is 
postponed to further work. 
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FIG. 5: (Color online) Emission pattern for the same configu- 
ration as in the right hand panel of Fig. 2] but in the presence 
of a disorder potential acting on the polaritons. The height 
of the disorder potential is in the 1 meV range. Inset: density 
profile of the disordered condensate. 

The case of a weak disorder is illustrated in Fig. [5l the 
strength of the disorder potential is chosen to be strong 
enough to substantially modulate the polariton density 
(see inset), but weak enough not to split the condensate 
into several frequency components [2^. In this case, the 
FWM signal remains clearly visible: the wave vector dis- 
tribution is strongly broadened by the disorder, but the 
main features that were observed in the right panel of 
Fig. [Hare still apparent. 



The situation is dramatically different whenever mul- 
tiple condensates are present with a substantial spatial 
overlap. In this case, the main effect of the probe beam is 
to redistribute the intensity among the different conden- 
sate frequencies and no clear FWM signal is easily identi- 
fied among the many spectral lines forming the emission 
spectrum. 



IV. CONCLUSIONS 

We have proposed and analyzed a four-wave mixing 
spectroscopy scheme to probe the excitation spectrum of 
a non-equilibrium polariton condensate. A clear evidence 
of the negative frequency 'ghost' branch is predicted to 
appear in the FWM spectra. The presence of this feature 
depends crucially on the coherent nature of the conden- 
sate and can therefore be used as an additional probe of 
the condensate coherence. Mapping out the FWM reso- 
nance as a function of frequency and wave vector gives a 
precise access to the elementary excitation spectrum. We 
have numerically demonstrated that a spatially localized 
probe can be used to overcome inhomogenous broaden- 
ing effects by restricting the measurement to the central 
region of the condensate where the density is almost flat. 
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